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Abstract 

We study the finite-size behavior of two-dimensional spin-glass models. We consider the ±J 
model for two different values of the probability of the antiferromagnetic bonds and the model 
with Gaussian distributed couplings. The analysis of renormalization-group invariant quantities, 
the overlap susceptibility and the two-point correlation function confirms that they belong to the 
same universality class. We analyze in detail the standard finite-size scaling limit in terms of TL l / u 
in the ±J model. We find that it holds asymptotically This result is consistent with the low- 
temperature crossover scenario in which the crossover temperature, which separates the universal 
high-temperature region from the discrete low-temperature regime, scales as T C (L) ~ L~ 6s with 
6 S « 0.5. 
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I. INTRODUCTION 

The two-dimensional (2D) Edwards-Anderson spin-glass model [1] has been extensively 
studied in recent years in order to investigate the interplay of disorder and frustration 
in 2D systems. If frustration is sufficiently large, these systems are paramagnetic at any 
finite temperature T. A critical glassy behavior is only observed for T — > 0. The zero- 
temperature behavior has been extensively studied. It is now well understood that it depends 
on the behavior of the low-energy spectrum. One should indeed distinguish systems with 
a discrete energy spectrum (DES), such as the ±J Ising model with a bimodal coupling 
distribution, from systems with a continuous energy spectrum (CES), such as the Ising 
model with Gaussian distributed couplings At T = 0, these two classes of systems 

behave quite differently. For instance, in DES systems, the stiffness exponent vanishes, while 
in CES systems, we have 9 < 0; recent numerical studies [5] give 9 ~ —0.28. 

For finite values of T, CES systems show a standard critical behavior, consistent with 
what is observed at T = 0. In particular, one expects v = — 1/9 and r\ = 0, two predictions 
which are consistent with numerical data 

a a 

The behavior of DES systems is instead 
more complex. In a finite box of linear size L, one observes two different behaviors, which 
depend on how large L is compared to a temperature-dependent crossover length L C (T); see 
Refs. B,Bfl and references therein. For L < L C (T), the critical behavior is analogous 
to that observed at T = 0. The system shows an effective long-range spin- glass order and 
its critical behavior can be predicted using droplet theory E|- On the other hand, for 
L > L C (T), the system is effectively paramagnetic. Equivalently, at fixed L, one observes the 
two regimes for T < T C (L) and T > T C (L), respectively, where T C (L) is the corresponding 
crossover temperature. Note that this discrete behavior can only be observed for finite 
values of L, since T C (L) — > for L — > oo. Of course, since T C (L) is an effective finite- 
size temperature, the crossover temperature is not uniquely defined, and many different 
definitions can be used. One of the basic questions is whether the critical behavior of 
DES systems for T > T C (L) is the same as that observed in CES systems. The numerical 
results of Ref. Q| strongly suggested that this is the case. However, those conclusions were 
later questioned in Ref. [6(, on the basis that much larger lattices were needed to show it 
conclusively. 

In a renormalization-group (RG) picture, the two regimes can be interpreted as due to 
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two different fixed points (FPs) a stable FP — the same that controls the critical 

behavior of CES systems — which describes the infinite- volume behavior up to T = 0, and 
an unstable FP, present only in DES systems, which controls the low-temperature behavior 
for T < T C (L). 

In order to fully specify the two regimes, one should predict how T C (L) scales with the 
size L. A free-energy argument, based on the energy difference and degeneracies of the 
two lowest-energy states p, suggests T C (L) ~ 1/ lnL as L — > oo. However, recently, using 
droplet theory, Refs. |7H9| suggested a power-law behavior 

L C (T)~T-^, T c (L)~L- e «. (1) 

Reference [3] predicted 6$ ~ 0.50(1), which appears to be consistent with their numerical 



8], as well as with 



data for the twist free energy [/J and the two-point correlation function 
previous results 4], 11]. A calculation in a hierarchical model js] gives a similar result of 
9s ~ 0.37. These calculations indicate that although 9s is quite small, it is nonetheless 
larger than the exponent 1/p — —9 « 0.28 (9 is the stiffness exponent in CES systems). 



In this paper, we investigate again the question of universality, by comparing the finite- 
size scaling (FSS) of the ± J model for two values of the disorder parameter p = 0.5, 0.8 and 
the model with Gaussian distributed couplings (henceforth we call it the Gaussian model). 
The FSS analysis in terms of RG invariant quantities (for example, the plots of the Binder 
cumulants versus the ratio with all quantities being defined in terms of the overlap 
variables) shows that the two models belong to the same universality class, confirming the 



conclusions of Refs. 



4|. Indeed, the ± J data have the same FSS behavior as the Gaussian 
data, if we only consider the ± J model results corresponding to temperatures larger than 
the crossover temperature. Then, we focus on the validity of the standard FSS in terms of 
the variable TL 1 ^, which is a rather subtle point in DES systems. Standard FSS exists only 
if T c (L)L l ' v -> for L -> 00. If we assume T C (L) ~ L 9s , then since T c (L)L 1//l/ ~ L 1 ^ 0S , 
FSS can be observed only if 6 S > 1/p w 0.28. This implies that if T C (L) ~ 1/ lnL flQ, the 
FSS limit T — > 0, L — > 00 at fixed TL 1 ^ does not exist in DES systems. On the other hand, 
if Eq. (CE]) holds with 9s ~ 0.50, FSS holds also in DES models. However, the approach to the 
asymptotic limit is quite slow. The region in which no FSS is observed, which corresponds 
to TL X I V < T C {V)L X I V , shrinks slowly, as L 1 ^" 03 ~ L~ a2 . The comparison of the Monte 
Carlo (MC) simulations of the ± J and Gaussian models shows quite convincingly that for 
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fixed TL X I V close to T C {V)L X I V ', the ± J model data converge toward the data of the Gaussian 
model, confirming the existence of the standard FSS, hence the power-law behavior (CQ) with 
9g > \jv. A reanalysis of the freezing temperature Tf(L), defined in Ref. lOj from the 
freezing of £/L and of the Binder cumulants [at fixed L, they are approximately constant 
for T < Tf(L)}, is consistent with Eq. ([[]) with 6$ ~ 0.4, which is close to the estimate 
of Ref. The freezing temperature T/(L) should represent a correct effective definition 
for T C (L), although deviations from the universal FSS behavior are expected for somewhat 
larger values of T. 

We also investigate the FSS behavior of the magnetization and the two-point correlation 
function of the overlap variables. We find that the data are consistent with the hyperscal- 
ing relation 2/3 = r\v. However, the data are not sufficiently precise to provide a precise 
determination of n, being consistent with a small value of ri < 0.2, including n = 0. In 
Ref. [8| , the authors showed that a properly subtracted overlap correlation function scales in 
the temperature region they consider, which essentially corresponds to T < T C (L). Here we 
consider the opposite regime, T > T C (L). We find that standard FSS as well as universality 
hold for the overlap correlation function. 

The paper is organized as follows. In Sec. HIl we define the models and the quantities 
we investigate. Section [TH] reports the numerical results of our FSS analysis: in Sec. IIII Al 
we discuss the RG invariant couplings, such as the ratio £/L and the cumulants of the 
overlap variable, focusing mainly on the question of the validity of FSS in terms of TL 1 ^; in 
Sec. IIIIBl we discuss the overlap magnetization and susceptibility; and finally, in Sec. IIII CI 
we discuss the two-point correlation function. In Sec. HVl we present our conclusions. 



II. MODELS AND DEFINITIONS 



We consider the 2D Ising model on a square lattice with Hamiltonian 

H= — / J Jxy&xCTy, (2) 

where a x = ±1, the sum is over all pairs of lattice nearest-neighbor sites, and the exchange 
interactions J xy are uncorrelated quenched random variables. We consider a model with 
Gaussian bond distribution, 

P(J xy ) ~ exp(- JIJ2) (3) 
4 



(in the following, we call it the Gaussian model). We also consider the ±J model where the 
couplings J xy take values ±J with probability distribution 

P(J xy ) = pS(J xy - J) + (1 - p)5(J xy + J). (4) 

As in Ref. [lo|], we consider p = 0.5 and 0.8. We recall that for sufficiently large frustration, 
i.e., for 0.11 < p < 0.89, the model shows a zero-temperature glassy critical behavior, 
with a paramagnetic low-temperature phase. Ferromagnetism can only be observed for 



p>p* = 0.89093(3) [12]. 



The critical modes at the glassy transition are those related to the overlap variable 
q x = a x lS) a x 2 \ where the spins ci belong to two independent replicas with the same dis- 
order realization {J xy }- In our Monte Carlo (MC) simulations, we compute the overlap 
magnetization 



4«l £&!>]' (5) 



m u 

X 

the overlap susceptibility x> an d the second-moment correlation length £ defined from the 
correlation function 

G o (x) = [(q q x )] = i(a a x )% (6) 

where the angular and the square brackets indicate the thermal average and the quenched 
average over disorder, respectively. We define x = J2 X G (x) and 

f2 _ 1 GM^-Gojp) 
4sin 2 (p min /2) G {p) 

where p = (p m i n , 0), p min = 2n/L, and G (q) is the Fourier transform of G Q (x). We also con- 
sider some quantities that are invariant under RG transformations in the critical limit, which 
we call phenomenological couplings. We consider the ratio £/L and the quartic cumulants 

TT - ^ TT - [Ml ~ If^ 2 (R\ 

where p k = ( Ej, ) k ). 

In the case of a T = transition with a nondegenerate ground state, as expected in 
CES systems, the overlap magnetization exponent vanishes, and C/4 — > 1 and U22 — > 
for T — > 0. Moreover, assuming the hyperscaling relation 2/3 = r]i>, we obtain rj = 0, thus 
X ~ £ 2 for T 0. 
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FIG. 1: (Color online) The cumulants U4 and U22 vs £/L. We present data for the ±J model 
(p = 0.5 and 0.8) and for the Gaussian model. Only the data for the largest lattices are included 
for clarity. 

III. FINITE-SIZE SCALING BEHAVIOR 



In order to study the FSS behavior, we extend the MC simulations of the ±J model 
at p — 0.5 and 0.8 presented in Ref. 10] ; we perform further simulations of the ±J Ising 
model at p — 0.8 on finite square lattices of sizes L = 16, 32, and of the Gaussian model for 
L = 8, 12, 16. We use the Metropolis algorithm and the random-exchange method [3]. For 
the Gaussian model, we study the temperature interval T min < T < 1.6, with T min = 0.2 
(L = 8) and T min = 0.167 (L = 12, 16). We average over a large number of disorder samples, 
i.e., 10 4 for each T and p. 



A. RG invariant couplings 



To begin with, we wish to check that the ±J model and the Gaussian model belong to 
the same glassy universality class, extending the FSS analyses of Refs. ay , Q. For this 
purpose, we consider [7 4 and U22 as a function of £/L. Our numerical results for the largest 
values of L are reported in Fig . [TJ No scaling corrections are visible in the plot of C/ 4 , as 
already observed in Refs. JfJ flQ| . while slightly larger corrections appear in the case of U22 for 
£/L > 0.6. It is, however, evident that as L increases, the differences between the Gaussian 
model results and those for the ± J model decrease. Thus, these results, together with those 
presented in Refs. {3, 6] (in Ref. js), other CES and DES systems were considered), confirm 
that all models belong to the same universality class. 
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FIG. 2: (Color online) The phenomenological couplings £/L (left) and U22 (right) vs TL X I V . We 
report the results for the Gaussian model (top), and for the ±J model at p = 0.5 (middle) and at 
0.8 (bottom). 

Now we investigate the question of the existence of standard FSS as a function of TL 1 ^ . 
For a RG-invariant quantity R, we expect 



R = hft(x 



x 



aTL l/l 



(9) 



where a is a nonuniversal constant that depends on the model but not on the quantity R, 
which can be chosen so that Hr(x) is model independent. 

First, we consider the data for the Gaussian model. The numerical results for £/L and 



U22 are reported versus TL X I V in Fig. [2] (top). We use v = 3.55, which corresponds to 
9 = —1/v = —0.282, which in turn is the present best estimate of the stiffness exponent 
of the Gaussian model jsj. The data show that £/L and U22 scale reasonably as a function 
of TL X I V and clearly appear to approach a FSS limit with increasing L. Scaling violations 
increase as TL 1 ^ increases. This is not unexpected since the correct scaling variable is 
the combination w<r(T)L 1 / !y , where ut{T) is the nonlinear scaling field associated with the 
temperature. Considering TL 1 ^ as the scaling variable corresponds to expanding ut(T) to 
first order in the temperature, an approximation which is expected to work well only when 
T is small. On the other hand, the region TL l / u > 2 corresponds to temperatures T > 1 
for the lattice sizes considered here. Of course, we cannot exclude the additional presence 
of nonanalytic scaling corrections, which increase as TL l / u increases. 

Then, we consider the data for the ± J model; see Fig. [2j The data show essentially 
three types of behavior, depending on the value of TL l / u '. For TL X I V < 1.5, no scaling is 
observed. This may be explained by the fact that the data in this region are below the 
crossover temperature, i.e., they correspond to T < T C (L), and thus are outside the regime 
in which FSS is supposed to hold. Then, there is an intermediate region, 1.5 < Tl}l v < 2.0, 
where data show scaling with small corrections — this is particularly evident for U22 ■ For 
TL X I V > 2.0, the corrections are larger, but the results appear to rapidly converge to a 
limiting curve: for both £/L and V221 the results satisfying L > 32 are very close to each 
other. We conclude that at least for TL 1 ^ > 1.5, FSS apparently holds. 

Furthermore, we verify the universality of the FSS behavior by comparing the results 
for the function Hr(x) defined in Eq. fl9]). For this purpose, we should first fix the model- 
dependent constant a that appears in Eq. (Q. We determine it by requiring the FSS curves 
for £/L to coincide in the region in which £/L ~ 0.5. Indeed, in this range of values of 
we observe small scaling deviations in all of the models we consider. If we set a = 1 for the 
Gaussian model, then for the ± J model, we obtain a(p = 0.5) ~ 1.3 and a(p = 0.8) ~ 1.5. In 
Fig. [3j we plot together the data for the Gaussian model and the ± J models at p = 0.8 and 
0.5. For clarity, we only report the data with L > 32 for the ± J model and the results with 
L > 12 for the Gaussian model. With this choice, there is only one point (it belongs to the 
± J model with p = 0.8 and corresponds to L = 32) which belongs to the region T < T C (L). 
This point is clearly visible in the figures as an isolated point. If we discard it, we observe 
good scaling up to aTL 1 ^ < 1.5: the ±J model data and the Gaussian data fall on top of 
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FIG. 3: (Color online) The phenomenological couplings £/L (left) and U22 (right) versus aTL 1 ^. 
We report the results for the ±J model at p = 0.8 and 0.5 (only L > 32) and for the Gaussian (G) 
model (only L > 12). We fix a = 1, 1.5, and 1.3 for the Gaussian model, and for the ±J model at 
p = 0.8 and 0.5, respectively. 

each other with good precision. For aTL 1 ^ > 1.5, the ±J data scale reasonably. The data 
of the Gaussian model, which correspond to significantly smaller lattices, show significant 
scaling corrections. It is, however, reassuring that the trend is correct: as L increases, they 
approach the ±J results. 

In order to understand the behavior close to the crossover temperature T C (L), in Fig. H] 
we report results for all values of L, but only for aTL 1 ^ < 1.25. It is clear that the 
deviations between ± J and Gaussian data slowly decrease as L increases (the same occurs 
for p = 0.5, not shown). This is consistent with the idea that FSS in terms of Tl}l v holds 
asymptotically, and, hence, with the prediction T C (L) ~ L~ 8s with 9s > 1/v. The approach 
is, however, very slow. Indeed, the region in which FSS does not hold is predicted to shrink 
as L- e s+V- „ L ^-\ 

This power-law behavior is also supported by the scaling of the freezing temperature 
Tf{L) defined in Ref. [h|. For each value of L, £/L and U4 become constant for small T, 
assuming values (£/L)f and U^j. Then, one defines Tf(L) as the largest temperature of the 



;o a power-law behavior 



151 ] used to determine Tf(L), it 



region in which £/L m (£/L)f and U4 ~ U±j. A fit of the date 
gives Tf(L) ~ L~ a35 ; see Fig. [51 Given the ad hoc procedure 
is difficult to give a reliable error for the result. It is, however, reassuring that the estimate 
satisfies the bound 6s > 1/u and is close to the estimates of Refs. 0, Q|. 
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FIG. 5: (Color online) Log-log plot of the freezing temperature Tt(L) vs L for the ±J model at 
p = 0.8. The dashed line shows a fit to the data, corresponding to Tf(L) ~ L~ 0,35 . 

B. Overlap magnetization and susceptibility 

The overlap magnetization m and susceptibility are expected to behave as 

m = C P/u u h (T)f m (Z/L), (10) 

and 



(11) 



where f m (x) and f x (x) are universal functions apart from a multiplicative constant, and the 
scaling field Uh{T) is an analytic function of T. If hyperscaling holds, we should have 



P 



(12) 
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FIG. 6: (Color online) Scaling combination H defined in Eq. (|13p vs £/L for the Gaussian (G) 
model and the ± J model at p = 0.8. 



Thus, the combination 



H 



£ 2 m 2 



(13) 



should be a universal function of independent of the scaling field Uh(T). In Fig. [6], 
we show the combination H for the Gaussian model and the ±J model at p = 0.8. The 
scaling is good. Deviations are only observed for £/L < 0.2 — these data correspond to 
large temperatures — and for £/L > 0.6, which, as discussed in Ref. lOi j . is the region 
in which strong crossover effects are observed for the lattice sizes considered in this paper. 
However, the observed trends are consistent with a unique universal curve. We thus confirm 
the validity of Eq. ( !T2|) . independently of what the numerical value of r\ is. If, indeed, r\ = 0, 
as theoretically predicted, Eq. f|T2|) gives (3 = 0. 

We have fitted all Gaussian data for the overlap susceptibility to Eq. (jl~T|) — more precisely 
to its logarithm as in Ref. [h3] — obtaining 77 = 0.20(7). The error we report is purely 
statistical and does not take into account possible scaling corrections. This result is slightly 
larger than the predicted result 77 = 0. The discrepancy should not be taken seriously, given 
the small lattices we consider. A precise determination of 77 in the Gaussian model requires, 
indeed, much larger values of L; see Ref. jfj]. 

The function f x (x) is universal apart from a rescaling: if fi(x) and f2{x) are determined 
in two different models, we expect f\(x) = bf'2 (x), where b is a model-dependent constant. 
We now compare the estimates of the functions f x {x) for the Gaussian model and the ±J 
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and 0.5. 

model. In Fig. [7J we report the functions (they have been rescaled so that they coincide for 
£/L — 0.3) for the Gaussian model and the ±J model at p — 0.8 and 0.5. We report the 
curves both for rj = and 0.2. In Ref. [10], we observed that for rj — 0, the scaling was good 
U P to £/L ~ ((/Ijmaxi a value which had been estimated as the boundary of the crossover 
region before the regime in which freezing was observed. The quantity (£/£) m ax should 
scale as (£/L) c , i.e., the value of £/L at the crossover temperature. Now, using Eq. (Q, 



since /i 



x) ~ x " for x — > to recover the correct infinite-volume behavior, we have 

(£/£) max ~ (£/L) c ~ ^[aT^L 1 /*] ~ T c (L)-"/L ~ L^- 1 ~ L - 8 , (14) 

where we have used = 0.5 in the last step, and the fact that T c (L)L l l u — )■ for L — )■ 00. 
Again note that the inequality > 1/z/ is necessary to guarantee that (£/L) ma x — ► 00 as 
L — y 00. By looking at the scaling behavior of t/22, Ref. [10] estimated (£/L) max ~ 0.65 
and 0.45 for p = 0.8 and 0.5, respectively, in the range 32 < L < 64. For the ±J model at 
p = 0.5, a similar estimate is obtained by considering the scaling behavior of the estimates of 
£(2L)/£(L) reported in Ref. [3(. The data for the Gaussian model agree with the ±J model 
data for both p = 0.5 and 0.8 up to (£/L) max . Thus, the numerical results are consistent 
with universality and rj = 0. 

In Ref. [lO |, we also observed that if we included all data in the fits, the best estimate of 
rj was 77 ~ 0.2. Indeed, in this case, the results for the overlap susceptibility showed a very 
good scaling up to £/L ~ 0.8. We did not take this result as an indication that 77 = 0.2 was a 
more plausible estimate than 77 = because we had good reasons to discard all data beyond 
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FIG. 8: (Color online) Scaling combination g Q (x,T) for x = L/2 (left) and x = L/4 (right) vs ^/L. 
We report results for the Gaussian (G) model and the ± J model at p = 0.8. 

£/L ~ (^/I/) max - Somewhat surprisingly, if we now include the data for the Gaussian model, 
and take 77 = 0.2 for all models, we again observe a good universal scaling up to £/L ~ 0.8. 



However, the numerical results of Ref. 



6| exclude rj = 0.2 for the Gaussian model. Thus, the 



apparently good observed behavior cannot hold asymptotically. In view of the result rj = 
of Ref. 6j, as L increases, the Gaussian data should become inconsistent with 7] = 0.2. 

Finally, we note that, in order to observe a good scaling behavior for the overlap sus- 
ceptibility, it is crucial to include the nonlinear scaling field Uh(T). Indeed, such a function 
gives a sizable contribution to our data. In the case in which we set rj = 0, we obtain 
u h (T = l) 2 /u h (T = 0.2) 2 ps 2.2 and u h {T = l.hf/u h (T = 0.2) 2 ps 3.3. 



C. Two-point function 

In Ref. 8|, the authors analyzed the scaling behavior of the two-point function, showing 
that below the crossover temperature, the two-point function scales as predicted by droplet 
theory. Here we analyze the two-point function in the opposite regime in which we expect 

G (r,T) = C 2 xf G (r/^Z/L). (15) 

Note that by writing the scaling behavior in this form, there is neither the need to specify 
r) nor to introduce the nonlinear scaling fields. Moreover, the function fc(x,y) is universal. 

To verify the scaling behavior (fl~5!) . we compute G D (r,T) along a lattice line, i.e. for 
r = (x,0). We perform the computation in the Gaussian model (L = 8, 12, 16) and in the 
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FIG. 9: (Color online) Scaling combination g Q (x,T) for £/L ~ 0.18 (left) and £/L ~ 0.63 (right) 
vs We report the results for the Gaussian (G) model and the ± J model at p = 0.8. 



±J model at p = 0.8 (L = 16, 32). Then, we consider 

G (r,T)L 2 



g a (x,T) 



\ 



(16) 



In Fig. El we report g {x, T) for x = L/2 and L/4 as a function of £/L, while in Fig. |9l we 
show g Q (x,T) at fixed £/L ~ 0.18 and £/L ~ 0.63 as a function of x. In these cases, the 
scaling is very good: all points fall onto a single curve, confirming the validity of Eq. ( IT5|) 
and universality. 

IV. CONCLUSIONS 

We investigate the FSS behavior of two-dimensional Ising spin-glass systems. In partic- 
ular, we consider the square lattice ± J model at p — 0.5 and 0.8, and the Gaussian model. 
In this respect, the ± J model appears particularly problematic because it presents two dif- 
ferent finite- volume regimes: a continuous regime for T > T C (L) and a discrete regime for 
T < T C (L). According to droplet theory, the crossover temperature T C (L) is expected to 
vanish in the large-L limit as a power law [l-^], T C (L) ~ L~° s with 6s ~ 0.5. A logarithmic 
behavior, T C (L) ~ 1/lnL, is instead suggested by the free-energy arguments of Ref. [6|. 

The main conclusions of our numerical analysis based on MC simulations are as follows: 

(i) All models we consider belong to the same universality class. The magnetization, sus- 
ceptibility, two-point correlation function, and the quartic cumulants, defined in terms 
of the overlap variables, show a universal FSS behavior in terms of £/L. In the case 
of the ±J model, this universal scaling is only observed above the crossover tempera- 
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ture T C (L), which separates the continuous region from the discrete low-temperature 
behavior. 

(ii) Our FSS analysis provides good evidence of the FSS limit T — > 0, L — > oo at fixed 
TL l l v . This implies that the crossover temperature T C (L) does not behave as 1/hiL 
as suggested in Ref. [6|, but rather as T C (L) ~ L~° s with 6 S > \jv « 0.28. This is 
consistent with droplet theory, which predicts a power-law behavior with 6s ~ 0.5. 

(iii) We study the FSS behavior of x- The data for the Gaussian and the ±J models 
support universality. However, the available data are not sufficient to obtain a precise 
estimate of rj and confirm definitely the expected value n = 0. 

(iv) We verify the hyperscaling relation 2/3 = r\v. If /3 = 0, it implies r] — 0. 

(v) We consider the two-point correlation function and show that it satisfies a stan- 
dard FSS ansatz. Note that the scaling form (fl~5l) . which is appropriate for a high- 
temperature phase, is different from that considered in Ref. which is appropriate 
for a low-temperature phase in which spontaneous magnetization is present. The oc- 
currence of these two different scaling behaviors is related to the different regimes 
considered. Here, we use data such that T > T C (L), while Ref. [8j studies the behavior 
for T < T C (L). 
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